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ABSTRACT 

This paper reports the discovery and characterization of the transiting hot giant exoplanet Kepler- 
17b. The planet has an orbital period of 1.486 days, and radial velocity measurements from the Hobby- 
Eberly Telescope show a Doppler signal of 419.5^![5g ms^^. From a transit-based estimate of the host 
star's mean density, combined with an estimate of the stellar effective temperature Teff=5630 ± 100 
from high-resolution spectra, we infer a stellar host mass of 1.06 ± 0.07 Mq and a stellar radius of 
1.02 ± 0.03 Rq. We estimate the planet mass and radius to be Mp — 2.45 ± 0.11 Mj and Rp = 
1.31 ± 0.02 i?j. The host star is active, with dark spots that are frequently occulted by the planet. 
The continuous monitoring of the star reveals a stellar rotation period of 11.89 days, 8 times the the 
planet's orbital period; this period ratio produces stroboscopic effects on the occulted starspots. The 
temporal pattern of these spot-crossing events shows that the planet's orbit is prograde and the star's 
obhquity is smaller than 15°. We detected planetary occultations of Kepler-17b with both the Kepler 
and Spitzer Space Telescopes. We use these observations to constrain the eccentricity, e, and find 
that it is consistent with a circular orbit (e < 0.011). The brightness temperatures of the planet the 
infrared bandpasses are T3.6 um=1880± 100 K and T4.5 ^ni=1770± 150 K. We measure the optical 
geometric albedo Ag in the Kepler bandpass and find Ag = 0.10 ± 0.02. The observations are best 
described by atmospheric models for which most of the incident energy is re-radiated away from the 
day side. 

Subject headings: planetary systems — stars: individual (Kepler-17b, KIC 10619192, 
2MASS 19533486-F4748540) — eclipses — techniques: photometry 



1. INTRODUCTION 



of Technology, Cambridge, MA 



^ Harvard-Smithsonian Center for Astrophysics, 60 Garden 
Street, Cambridge, MA 02138; jdesert@cfa.harvard.edu 

^ Massachusetts Institute 
02159, USA 

^ Department of Astronomy and Astrophysics, University of 
Cahfornia, Santa Cruz, CA 95064, USA 

* University of Texas, Austin 

^ Department of Astronomy, University of California, Berke- 
ley CA 94720-3411, USA 

° Solar System Exploration Division, NASA Goddard Space 
Flight Center, Greenbelt, MD 20771, USA 

^Neils Bohr Institute, Copenhagen University 

* NASA Ames Research Center, Moffett Field, CA 94035 
® San Jose State University 

^^ Las Cumbres Observatory Global Telescope, Goleta, CA 
93117 

" SETI Institute/NASA Ames Research Center, Moffett 
Field, CA 94035 

1^ University of Florida, Gainesville, FL 32611 

^^ Space Telescope Science Institute, Baltimore, MD 21218 

^"^ Institut d'Astrophysique et de Geophysique, Universite de 
Liege, AUee du 6 Aofit 17, Bat. B5C, 4000 Liege, Belgium 

Bay Area Environmental Research Inst. /NASA Ames Re- 
search Center, Moffett Field, CA 94035 

^^ Orbital Sciences Corporation/NASA Ames Research Cen- 
ter, Moffett Field, CA 94035 

^'^ Centre for Astrophyiscs Research University of Hertford- 
shire 



NASA's Kepler mission is a space-based photomet- 
ric telescope dedicated to finding Earth-size planets 
in the habitable zones of their host stars and to de- 
termining the frequency and char acteristics of plane - 
tary systems around Sun- like stars ( Borucki et al.l l2010) . 
While monitoring nearly continuously the brightness of 
about 150,000 dwarf stars since May 2009 to achieve 
the above goals, it accumulates an extraordinary volume 
of data allowing atmospheric studies of the transiting 
hot-J upiters present in the field of view (jBorucki et al.l 
l2009f l. The first four months of observations were re- 
leased and 1235 transiting planet candidates were re- 
ported (jBorucki et al.|[20Tll ). Because of their short pe- 
riods and relatively deep transit depths, several hot- 
Jupiter candidates were discovered early on during the 
first weeks of the mission. Amongst those, the Kepler 
Object of Interest KOI-203 was identified as a promis- 
ing target, was selected for follow-up studies, and placed 
on the short-cadence target list for subsequent quarters. 
This is the object of interest of the current study. 

The highly irradiated transiting hot-Jupiters currently 
provide the best opportunities for studying exoplane- 
tary atmospheres in emission, during planetary occulta- 
tions, when the exqplanets pass behind th e ir parent stars 
(iSeager et all 120001: iSudarskv et all 120001: HbTtnev et al.l 
120051 : iBarman et al.l 120051^ Light emitted from exo- 
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planets was first detected froni space at infrared wave - 
lengths (jCharbonneau et all 120051 : Deming et al. 1200 " 



and more recently i n the optical (lAlonso et al 



2009 



Snellen et all 120091 : iBorucki et al.l 120091 : lAlonso et all 



2010( ) us ing the CoR oT an d Kep ler Space Telescopes. 
Notably, iRowe et al.l (|2006L 120081 ) used the Microvari- 
ablity and Oscillations of Stars {MOST) telescope to 
place a very stringent upper limit on the depth of the 
occultation of HD 209458b. 

The hot-Jupiters detected by CoRoT and Kepler are 
particularly good targets for studying planetary atmo- 
spheres, because the precise and nearly uninterrupted 
photometric surveillance provided by the space satellites 
can allow the planet ary occultatio i is to be detected at 
optic a l wavelengths (iSnellen et al.l 120091 : IBorucki et al.l 
120091 : iDesert et all l2011bD . Obtaining multiple wave- 
length observations of the relative depths of planetary 
occultations is necessary to constrain the broad band 
emergent spectra. Such observations are fundamen- 
tal to underst anding the energy budget of the s e ob- 
jects ( Sudars kv et al.l [2Q0l iBurrows et all [20051 [20081: 
iSpiegel k. Burrowsl 120101 ) and for comparative exoplan- 
etology. 

We confirm here the planetary status of KOI-203, des- 
ignating it as Kepler-17b, and study its atmosphere. We 
first describe the Kepler observations and transit model- 
ing as well as the follow-up observations used to confirm 
this planet, including an orbital solution using radial ve- 
locities obtained with the High Resolution Spectrometer 
(HRS) at the Hobby-Eberly Telescope (HET). We then 
study the impact of the stellar variability on the transit 
lightcurves and we use occulted dark starspots to con- 
strain the stellar obliquity. We finally combine occulta- 
tion measurements obtained in the optical with Kepler 
and in the infrared with Spitzer to learn about the at- 
mospheric properties of Kepler-17b. 

We first describe the observations, time series and anal- 
ysis of the Kepler photometry in Sect. 2, then we describe 
the follow-up observations that confirm the planet detec- 
tion in Sect. 3 and then discuss the stellar properties in 
Sect. 4. In Sect. 5, we describe the occultation measure- 
ments obtained from the visible and infrared. We present 
a global Monte- Carlo analysis of the complete sample of 
observations in Sect. 6 and finally discuss our findings in 
Sect. 7. 

2. Kepler OBSERVATIONS 

Observations of the Kepler field commenced in May 
2009 with Quarter (QO); the data that we describe here 
are the Kepler science data of Kepler- 17 from Quarter 
to 6 (Q0-Q6). The Kepler observations were gathered 
almost continuously during 16.7 months. These obser- 
vations have been reduced and detrended by the Ke- 
pler pipeline ([Jenkins et al.ll20l"ol ) . The Kepler bandpass 
spans 42 3 to 897 nm for which the response is greater 
than 5% (jBatalha et al.ll2010t [B7vson et al.ll2010[) . This 
wavel ength domain is roughly equivalent to the V-l-R- 
band (iKoch et al.ll2010[ ). The target Kepler- 17 was iden- 
tified in the Kepler Input Catalog (jBrown et al.l 120111 ) 
(2MASS 19533486-K4748540, KIC 10619192, r= 14.08 
mag). Because of the transit like events, the object 
was then considered for follow-up studies, identified as 
Kepler Object of Interest KOI-203. The pipeline pro- 
duces both calibrated light curves (PA data) for individ- 



ual analysis and corrected light curves (PDC) which are 
used to search for transits. This paper presents results 
which are measured from PA data only. They consist of 
long cadence i ntegration time (29.4 2 6 minutes) for Quar- 
ters and 1 fCald well et"all 120101 : [Jenkins et all |2010() 
an d long and short cad ence (1 minute) for Quarters 2 to 
6 ([Gilliland et al.|[2010[ ). The pipeline provides time se- 
ries with times in barycentric corrected Julian days, and 
flux in photo-electrons/per cadence. The raw nearly con- 
tinuous photometry of Kepler-17 is presented in Figure[TJ 
We measure the transit parameters from the Kepler ob- 
servations as described below (see Section[6|). We present 
the normalized, phase-folded and combined transit light 
curve obtained at short cadence in Figure [2] from which 
we measure the transit parameters. 

3. FOLLOW-UP OBSERVATIONS 
3.1. Reconnaissance spectroscopy 

As described in detail by iGautier et al.l ([2010I ). the 
follow-up observations of Kepler planet candidates in- 
volve reconnaissance spectroscopy to look for evidence of 
astrophysical false positives responsible for the observed 
transits. These false positives include single- and double- 
lined binaries, certain types of hierarchical triples and 
even some background eclipsing binaries, which would 
show velocity variations and/or composite spectra that 
are readily detectable by the modest facilities used for 
these reconnaissance observations. As described below, 
we also use these spectra to estimate the effective tem- 
perature, surface gravity, metallicity, and rotational and 
radial velocities of the host star. 

On 25 Apr 2011 UT, we obtained a spectrum of Kepler- 
17 using the fiber- fed Tillinqhast Reflector Echelle Spec- 
trograph (TRES; Fiircsz (2008)) on the 1.5m Tilling- 
hast Reflector at the Fred Lawrence Whipple Observa- 
tory., on Mt. Hopkins, AZ. The spectrum was taken 
with the medium fiber, which has a resolving power of 
A/AA=44,000 and a wavelength coverage of about 3850 
to 9100 Angstroms. The exposure time was 80 minutes. 
The spectrum was extracted and rectified to intensity 
vs. wavelength using standard procedures described by 
iBuchhave et al.l ([20101 ). The extracted spectrum has a 
signal-to-noise ratio of 17.5 per resolution element. 

We performed cross-correlations against a grid of syn- 
thetic stellar spectra, calculated by John Laird for a grid 
of Kurucz model atmospheres ([Kuruczl I1979D , using a 
line list developed by Jon Morse. The grid is coarse 
- 250 K in effective temperature, Tes; 0.5 in log sur- 
face gravity, log 5; 0.5 in log of the metallicity compared 
to the sun, [M/H]; and 2 km s-1 in rotational velocity, 
Vrot ~ so rather than simply adopting the parameters 
from the template with the best correlation coefficient, 
we fit a surface to the correlation peak heights to ar- 
rive at a refined classification. However, given the de- 
generacies between Tes, log 5, and [M/H], the quality 
of the spectrum is not sufficient to determine all three 
at once. As such, we have fixed \ogg — 4.42 ± 0.015 
as inferred from the best-fit of the Kepler photometry. 
The analysis yields roff=5630± 100 K, [M/H]=0.3±0.i, 
usini = 4.7±1.0 kms~^. When corrected for the orbital 
motion of Kepler-17 and the TRES zero point offset, de- 
termined by long-term monitoring of the lAU RV stan- 
dard HD 182488, we find the absolute mean systemic ve- 
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locityofKepler-17tobe-23.82±0.10kms"^ Note that 
this does not include any uncertainty in the absolute ve- 
locity of HD 18248 8, which we take t o be -21.508 kms"!, 
as observed by Nid ever et al.l (|2002| ). 

We also obtained Keck HIRES spectra and esti- 
mated the line strengths of 5hk=0.322 ± 0.01 and 
logi?^j^=-4.61 for Kepler-17 (assuming B-V=0.82). 
The Ca II H & K li ne strengths are a good indicator 
of the stellar activity (jlsaacson fc: Fischeiil2010[) . 



3.2. Imaging 

The astrometry derived from the Kepler images them- 
selves, when combined with high-resolution images of 
the target neighborhood, provides a very powerful tool 
for identifying background eclipsing bi naries blended 
with and contaminating the target images (jBatalha et al.l 
120101) . The astrometry of Kepler-17 indicated no signifi- 
cant offset during transits in any quarter, and computed 
offsets are well within the formal 3-sigma radius of confu- 
sion. Therefore Kepler-17 is considered to be the source 
for the transits observed in the Kepler lightcurves. 

We obtained an I-band image at the Lick Observatory 
1-meter, Nickel telescope with the Direct Imaging Cam- 
era (see Figure 131). The 1.0" arc-second seeing revealed no 
companions from 2" to 5" from the star's center, down 
to a limit of 19th magnitude. Similar conclusion were 
reached using UKIRT J-band images (Figure |3l). 

3.3. Radial velocity 

We obtained precise radial velocity (RV) follow-up 
observations of Kepler-17 with the HET (Ra msey et 
al. 1998) and its HRS spectrograph ()Tulllll998[ ) at Mc- 
Donald Observatory. Kepler-17 was observed ten times 
in the 2010 observing season, from 2010 August 22 un- 
til 2010 November 22. The instrumental set up and ob- 
servin g mode are described in more detail in lEndl et al.l 
(|201lD . Kepler- 17b is t he second planet confirmed with 
HET after Kepler-15b (lEndl et al.ll2"oTll) . We employed 
a "snap shot" strategy, using relatively short exposures 
of 1200 seconds, that yield a SNR sufficient to detect the 
radial-velocity signal of a hot- Jupiter. Thirteen spectra 
were taken with the l2-cell in the light path to compute 
precise differential RVs. These spectra have a typical 
S/N-ratio of 32 per resolution element. The radial veloc- 
ity data are listed in Table [TJ 

We use Gaussfi t , the g eneralized least-squares software 
of IJefFervs et al.l ((1988) to fit a Keplerian orbit to the 
HRS radial velocity data. Only the velocity zero-point 
and the radial velocity semi-amplitude K are included 
as free parameters in the fitting process. We first fitted 
the radial-velocity data alone, requiring the orbit to be 
circular (e = 0) and adopting the ephemeris derived from 
the Kepler photometry. The best-fit orbit has a K oi 
420 ± 15ms~^, a Xrcd ^^ 0-^ and a residual rms scatter 
around the fit of 52 ms~^. The radial velocity data and 
the orbital solution are shown in Figure |4l 

We determined the spectral line bisectors, which are 
a measure of line asymmetry, from the HET spectra to 
test if the radial velocity variations could be caused by 
distortions in the spectral line profiles due to contam- 
ination from a nearby unresolved eclipsing binary. We 
can only use a small fraction of the available spectral 
range that is not contaminated by the iodine absorp- 



tion cell (5000-6400 Angstrom) and thus the uncertain- 
ties in the bisector velocity span (BVS) are quite large 
with an average uncertainty of 99 ms^^. The RMS of 
the bisector measurements is 146 ms^^. There is no 
evidence of a correlation between the velocities and the 
bisectors, which supports the interpretation that the ve- 
locity variations ar e due to a planetary companion (e.g. 
lOueioz et aTllMll ) 

4. STELLAR PARAMETERS 

We derive the mass, radius, and age of the hos t star 
using the method described bv lTorres et al.l ()2008D . We 
first created a set of stell ar evolu tion models from the 
Yonsei-Yal e (Y^) series bvIYi et a l. (2001), with correc- 
tions from IDemargue et al.l (j2004f ). We employed their 
interpolation software which accepts as inputs the age 
of the star, the iron abundance, and the abundance of 
a-elements (relative to solar) for which we assume the 
solar value, and outputs a grid of stellar isochrones cor- 
responding to a range of masses. We evaluated a set of 
isochrones at ages from 0.1 to 14 Gyr (at intervals of 0.1 
Gyr) and stellar metallicities spanning a range of 3ct (at 
intervals of 0.01 dex) from the best-fit metallicity derived 
from spectra of [Fe/H]=0.3 ± 0.1. We then performed a 
spline interpolation of each output table at a resolution 
of 0.005 in effective temperature T^s, the log of the sur- 
face gravity \og{g), and the stellar luminosity L^. We 
evaluated the physical radius corresponding to each stel- 
lar model via log{g) and the mass of the star, though it 
is also possible to convert to physical radius using the 
model stellar luminosity and effective temperature (as- 
suming Li, = ATrR^aT^); in practice these conversions 
give identical results. 

We fitted for the stellar mass and radius us- 
ing Newton's versio n of K epler's third law in th e 
manner emplo yed by iSeaeer fc MaUcn-Ornelai (|2003| ). 
iSozzetti et all (|2007() and iTorres et al.l ([2008). We as- 
sumed that the planetary mass is negligible when com- 
pared to the mass of the host star. Using the Markov 
Chain Monte Carlo sequence of a/Ri,, and generating a 
series of Gaussian random realizations of [Fe/H] and Tog 
using the values and error bars derived from spectroscopy 
of Toff=5630± 100 K and [Fe/H]=0.3± 0.1, respectively, 
we located the best isochrone fit at each realization using 
the x^ goodness-of-fit given in Equation [TJ 



2 /Aa/i?,\' , /^ATeff' ' 
X = 



^a/R, 



V ^Teff 



A [Fe/H] 
.^[Fe/H], 



(1) 



Using the output of the MCMC chain of a/ Ri, ensures 
that any correlations between parameters, which are re- 
fiected in the chain, are properly incorporated into our 
estimate of the stellar parameters. 

We then assign a weight to the likelihood of each stellar 
model in the chain, applying a prior for the initial mass 
function (IMF) which assumes a Salpeter index. The 
number of stars of each mass and age, per 1000 stars, 
is generated by the interpolation software provided by 
lYi et al.l ()2001[ ) for several IMF assumptions, including 
the Salpeter IMF. We designate the weight assigned to 
each stellar model in the chain by normalizing to the 
highest IMF value within the sample; in practice, the 
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weights vary from 0.3 to 1 (from the least to most likely). 
The IMF prior changes the final answer by less than half 
a a for all parameters. We then incorporate this likeli- 
hood by discarding members of the chain according to 
their weight, where the weight is equal to the likelihood 
of remaining sampling; about 40% of the original chain 
remains intact after this stage. The value for each stel- 
lar parameter is then assigned from the median of this 
weighted distribution, with the error bars assigned from 
the nearest ±34% of values. In this way, we find Mj,= 
1.06ltEJ|^4g Mq, R^= 1.019 ± 0.014 Rq, and an age = 
2.9titg Gyr. These uncertainties are statistical as they 
exclude possible systematic uncertainties in the stellar 
models. The best-fit solution is presented in Figure [S] 

We caution that stellar isochrones are only poorly con- 
strained for this faint star. The estimated error bars on 
the stellar parameters are smaller than what is expected 
for such a star of this magnitude. This is because we 
set the stellar gravity as a fixed value in our analysis, 
since the stellar spectrum has a low signal-to-noise. We 
set the stellar gravity to the value we derive from the 
stellar density measured using the Kepler photometry 
(see Section [3T|) . The quality of the combined transit 
lightcurve allows us to measure the stellar density with 
a high precision that tightly constrains the isochrone fits 
(see Figure [S]). However, a recent study of the physical 
properties o f the stellar compon ents of the transiting ex- 
oplanets by iSouthwortbl ()2011[ ) shows that the current 
stellar models are not determined at better than 1% in 
terms of radius and 2% in mass (median values); more 
conservatives values are 3 and 5% for the stellar radius 
and mass respectively. We adopt these conservatives val- 
ues and propagate the systematic errors arising from the 
dependence on stellar theory to the final errors. Adding 
the statistical and systematical error quadratically, we 
find M^= 1.061 ±0.07 Mq, R^= 1.019 ± 0.03 Rq. This 
study shows that Kepler- 17 is consistent with a main- 
sequence G dwarf star, its radius and T^s are indistin- 
guishable from Solar. 

5. OCCULTATIONS OF Keplcr-17b 
5.1. Occultations from the Kepler photometry 

The first method we use to measure the occulta- 
tions from the Kepler observations is described by 
[Desert et al.l (|2011b( ). We search for the occultation 
events in the short cadence light curves. We normal- 
ize and combine 173 occultations events. The best-fit 
and the maximum depth as a function of the orbital 
phase are found very close to the orbital phase of 0.5 
as expected for a circular orbit. We estimate the sig- 
nificance of this detection by measuring the occultation 
depth, ephemeris and associated errors using a bootstrap 
Monte-Carlo analysis. We find that the planet has an oc- 
cultation depth of 58 ± 10 ppm. 

In a second method, we include the occultations in a 
Markov Chain Monte Carlo global analysis as described 
below in the Section |6l We present the normalized, 
folded, combined and binned per 15 min light curve ob- 
tained at long cadence in Figure[B]from which we measure 
the occultation depth with the MCMC. 

5.2. Spitzer observations and photometry 



The method we use to measure the occu hations from 
the S pitzer observations is described in iDesert et al.l 
(|2009D . Kepler- 17 was observed during four occultations 
between August 2010 and Januar y 2011 with Warm- 
Spitzer /IRAC (iWerner et al. l'2004':'F a^io et al.l[200l at 
3.6 and 4.5 fim (program ID 60028). Two occultations 
were gathered per bandpass and each visit lasted approx- 
imately 8.6 h. The data were obtained in full-frame mode 
(256 X 256 pixels) with an exposure time of 30.0 s per im- 
age and yielded 2461 images per visit. The method we 
use to produce photometric time se ries from the images 
is described in IDesert et al.l (|2011a[ ) . It consists of find- 
ing the centroid position of the stellar point spread func- 
tion (PSF) and performing aperture photometry using 
a circular aperture on individual Basic Calibrated Data 
(BCD) images delivered by the Spitzer archive. The op- 
timal apertures were found to be 2.5 pixels in radius. 
The final photometric measurements used are presented 
in Table [5] The raw time series are presented in the 
top panels of Figure [T] Telescope pointing drift results 
in fluctuations of the stellar centroid position, which, in 
combination with intra-pixel sensitivity variations, pro- 
duces systematic noise in the raw light curves (upper 
panel Figure [7]). A description of this effect, known 
as the pixel-pha se effect, is given i n the Spitzer /IRAC 
data handbook ([Reach et al.l l2006f) and is well known 
in exoplanetary stu dies (e.g. iCharbonneau et all 120051 : 
iKnutson et al .1120081 ). To correct the light curve, we de- 
fine a baseline function that is the sum of a linear func- 
tion of time and a quadratic function (with four parame- 
ters) of the X and Y centroid positions. We find that the 
point-to-point scatter in the photometry gives a typical 
signal-to-noise ratio (S/N) 90 and 110 per image at 3.6 
and 4.5 /im respectively. These correspond to 85% of the 
theoretical signal-to-noise. 

We simultaneously fit the instrumental functions with 
all the parameters, measured the occultation depths for 
each individual visit, and report report the values in Ta- 
ble [2] The measurements per bandpass agree at the 1- 
a level. The weighted mean of the transit depths are 
0.250 ± 0.030% and 0.310 ± 0.035% at 3.6 and 4.5 fim 
respectively. 

We also included the Spitzer observations in the global 
MCMC analysis described below. 

6. MARKOV CHAIN MONTE CARLO ANALYSIS 

The analysis of the Kepler photometry and the de- 
termination of the stellar and planetary parameters for 
Kepler- 17 follows pro cedures similar to those reported in 
iBorucki et al.l (|2010|) . We check and confirm these re- 
sults by an independent analysis with the global Markov 
Chain Monte-Carlo (M CMC) algorithm presented in 
iGillon et al.l ([20091 . [20To[ ) . This MCMC implementation 
uses the Metropolis-Hasting algorithm to perform the 
sampling. Input data to the MCMC include the Kepler 
and Spitzer photometric measurements and the radial- 
velocity data. 

Kepler Short Cadence (SC) data allow a precise deter- 
mination of the transit parameters and more specifically 
allow a fit for the limb-darkening (LD) coefficients. We 
therefore assumed a quadratic law and used ci = 2ui+U2 
and C2 = ui — 2u2 as jump parameters, where ui and U2 
are the quadratic coefficients. Those linear combinations 
help in minimizing correlations on the uncertainties of ui 
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and U2 (|Holmanet al.ll2OO60 . 

The MCMC has the following set of 13 jump pa- 
rameters: the planet-to-star flux ratios in the Kepler 
and Spitzer bandpasses, the transit depth, the impact 
parameter b, the transit duration from first to fourth 
contact, the time of minimum light, the orbital pe- 
riod, K' = Ky/1 — e^P^/"^, where K is the radial- 
velocity semi-amplitude, the two LD combinations Ci 
and C2 and the two p arameters ^/e cos uj and ^/esinu! 
([Anderson et al.ll2011f) . A uniform prior distribution is 
assumed for all jump parameters. Baseline model coeffi- 
cients ar e determined for e ach lightcurve with the SVD 
method (|Press et all 119921 ) at each step of the MCMC 
Corre lated noise is accounted for following I Winn et"aLl 
(|2008f ) and 'Gillon et all (|2010[ ). to ensure reliable error 
bars on the fitted parameters. For this purpose, we com- 
pute a scaling factor based on the standard deviation 
of the binned residuals for each lightcurve with different 
time bins. The error bars are then multiplied by this 
scaling factor. We obtained a mean scaling factor of 1.19 
for all photometry. The results from the MCMC analysis 
are presented in Table [3] 

7. RESULTS AND DISCUSSION 
7.1. Characteristics of the system Kepler- 17 

Adopting the stellar mass and radius for Kepler- 17 
that we found in Section |4l we obtain a mass of Mp = 
2.45 ± 0.11 Afj and a radius of Rp = 1.31 ± 0.02 Rj for 
the planet, which leads to a density of pp — 1.35 ± 
0.08 gcm~^, and a surface gravity loggp = 3.55 ± 0.02 
(cgs). We note that the surface gravity can be de- 
rived from the photom etry and the velocimetry only 
(jSouthworth et al.|[2007l) . The position of Kepler- 17b on 
the mass/radius diagram is fairly common and it appears 
to be slightly infiated compared to models that includ e 
the effects of stellar irradiation (e.g. iLatham et al.l[2010l ). 

The current upper limit on the orbital eccentricity e 
from radial velocity measurements is consistent with zero 
for Kepler-17b. We also measure the mid-occultation 
timing offset from both Kepler and Spitzer observations. 
The determination of the timing of the secondary eclipse 
constrains the planet's orbital eccentricity. Our estimate 
for the best-fit timing offset translates to a constraint on 
e and the argument of pericenter ui. The timing is used 
to constrain the ecos(a;) and we find that e < 0.011 at 
the 1 — a level. This upper limit implies that the orbits 
of these objects are nearly circular unless the line of sight 
is aligned with the planet's major orbital axis, i.e. the 
argument of periapse ui is close to 90° or 270°. 

We fit a linear function of transit epoch and find 

rc(0) = 2455185.678035i^ ;;g|]Ej26 (BJD„fc) with a period 
P = 1.4857108 ± 0.0000002 days. We also fit each mid- 
transit time individually and compare each one to the 
expected linear ephemeris to obtain a Observed-minus- 
Computed (0-C) time series. The 0-C times shows a 
scatter that is most probably caused by stellar spot- 
induced shifts, since the star is active (see Section [7^ . 
Formally, the individual fits are consistent with the lin- 
ear model, therefore we do not consider this to be a sig- 
nificant detection of timing anomalies. Since there is 
no clear evidence for transit timing varia tions (TTVs; 
lAgol et all [20051: iHolman fc Murravi[2005[ ). we use the 
timing data to place upper limits on the mass of a hy- 



pothetical second planet that would perturb the orbit of 
the transitiiig plan et using the procedure described in 
iCarter et al.l ([20 lit ). Figure [8| shows the constraints on 
the perturber mass as a function of period ratio, as de- 
termined from this analysis. The mass constraints on the 
perturber are more restrictive near the mean-motion res- 
onances and most restrictive at the low-order resonances, 
particularly for the interior and exterior 2:1 resonances. 
For example, a perturber at the interior 2:1 resonance 
having a mass near that of Mars would have induced 
TTVs detectable in the present data. 

7.2. Stellar variability of Kepler- 17 

The Kepler photometry exhibits a quasi periodic fiux 
modulation of about 3% with a period of 11.9 days (Fig- 
ure [1]). This period corresponds to 8 times the planet's 
orbital period (Figure [Hj) . The stellar variability can be 
quantified from a Lomb-Scargle periodogram that reveals 
two peaks at 5.95 and 11.9 days (see Figure ITUl) . The 
position of the peak on the periodogram and its width 
provide a measurement of the stellar rotation period and 
its error, respectively. We measured a stelUar rotation 
period of 11.9 ± 1.1 days. The presence of active stellar 
spots localized on opposite hemispheres can best explain 
the two peaks seen in the periodogram. We find that the 
ratio between the planet's orbital period and the stel- 
lar rotation period is 8.0 ± 0.7 days. We note that this 
intriguing integer ratio of 8 could potentially reveal the 
signature of stellar-planet interactions. 

7.2.1. Using occulted starspots to infer the orbital obliquity 

The lightcurve of Kepler- 17 shows substantial deforma- 
tion of the planet's transit profiles which we interpret as 
occultations of dark starspots (Figure [5]). We co-added 
the transit lightcurves that occurred at epochs modulo 
the stellar rotation period (modulo 8 planetary orbital 
periods). We present the resulting 8 transit lightcurves 
in Figure 1111 each contains about 22 individual transit 
lightcurves. We identified that the planet crosses spots 
of identical shapes every 8 transits producing a strobo- 
scopic effect. Therefore, we conclude that starspots come 
back to the same position on the transit chord after 8 or- 
bits or 11.89 days (one stellar rotation period). Because 
of this stroboscopic effect, the scatter measured during 
the spot-crossing phase of each of these lightcurve resid- 
uals is similar to the scatter of the residuals measured 
outside the transit. This further supports the assertion 
that the same spots are crossed by the planet every 8 
transits, at identical orbital phases, because the residu- 
als measured from a transit light curve obtained from a 
planet occulting a random distribution of spots are ex- 
pected to be larger during the transit event compared 
to the out-of-transit monitoring. This can be explained 
only if spots are coming back to similar longitudes and 
latitudes after one stellar rotation period. Interestingly, 
occulted stellar spot s can place constraints on the spin- 
orbit alignment (c.g Deming et al. 2011; Nutzman et al 
2011: Sanchis-Ojeda et al. 2011a: Sanchis-Oicda &: Winn 



120 im 

We identified five spots (A, B, C, D and E) by their 
anomaly in the residuals and we marked their initial or- 
bital phase positions as seen in Figure [11] We then 
computed the expected phase position of each of these 
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5 anomalies on the 8 following lightcurve residuals as- 
suming that occulted starspots should have moved by 
360/8 = 45° in longitude towards the egress (for a pro- 
grade and aligned orbit) for every planetary transit. We 
find that anomalies are always present at all the expected 
phases. 

As the star rotates, the spot surfaces projected on 
the stellar sphere change; therefore the shapes and am- 
plitudes of these anomalies change between consecutive 
transits while the same spots are crossed over time. Since 
a spot surface appears larger at the center of the star 
compared to the limb, the anomalies due to spots should 
exhibit larger amplitudes when they are crossed close 
to the mid-transit time. This is consistent with what 
is observed in Kepler- 17. For example, anomaly 'B' in 
Figure [Tl] is observed in four consecutive transits. The 
overall shape of this anolmaly is conserved transit after 
transit and its amplitude increases before and decreases 
after the mid-transit. 

We conclude that the planet crosses the same spots, 
transit after transit, at longitudinal positions that differ 
by 45 °. 

These findings imply that the projected spin-orbit an- 
gle, A, is very close to for this system. Because occulted 
spots can be identified transit after transit (Figure [T^ . 
with a change in phase expected from the stellar rotation 
period, we conclude that the stellar inclination angle is 
near 90 degrees. This is further supported by the mea- 
sured wsini = 4.7± 1.0 km/s which is in good agreement 
with its expected value of 2T:Ri,/Prot = 4.3 km/s assum- 
ing an inclination angle of 90 degrees. This implies that 
the true obliquity, i.e. the angle between the stellar ro- 
tational axis and a line perpendicular to planet's orbital 
plane, ^, must be close to 0. Finally, we conclude that 
the planet orbit is prograde since the occulted starspots 
progress from the transit ingress to the egress limbs. 

Assuming that the uncertainty on latitude of a spot 
corresponds to the planet size, we derive an uncertainty 
for * of ±arctan(i?p/i?*) = ±10°. We note that the 
spot size could be larger than the planet size, but this 
changes the results only modestly. Assuming a star spot 
twice the projected size of the planet on the stellar photo- 
sphere provides a more conservative uncertainty of ±15° 
fo r 'F. 

I Winn et al.l ()2010f) show that stars with transiting 
planets for which the stellar obliquity is large are prefer- 
entially hot (Teff > 6250 K). The low'obhquity of Kepler- 
17b fits this pattern since the stellar temperature is lower 
than 6250 K. 

7.2.2. Starspots Ufettme 

To help visualize the spot crossings during transit, we 
subtract the best-fit model from each individual transit 
and investigate the residuals. We slide a box of dura- 
tion twice the ingress time across the residuals in each 
transit epoch, recording the scatter for a given box mid- 
time (relative to mid-transit) and a given epoch. We 
then produce an image of this scatter at orbital phase 
as a function of the epoch number (see Figure [T2|) . We 
used an interpolation (with IDL's TRIGRID) with an 
output sampling of 600 samples in the epochal direction 
and 300 samples in the transit phase direction. This 
image reveals individual spots that we define as either 
"hot" or "cold" regions, depending on whether the indi- 



vidual slide box residuals are below or above the tran- 
sit light curve model. The repeating vertical structure 
is interpreted as spots marching across the transit chord 
such as seen in the previous section. Each vertical profile 
is slanted slightly from left to right indicating that the 
spots progress from the ingress limb to the egress limb. 
Some spots make their way around the star and reappear 
again during several stellar rotation periods. For exam- 
ple, the collection of "cold" spots in the image starting 
around Epoch 110 and ending around Epoch 170 seem 
to be related to the same spot. We conclude from the 
nearly continuous monitoring of Kepler- 17 that the oc- 
culted starspots are present on the same stellar chord for 
at least 100 days, somewhat comparable to the lifetime 
of sunspots. 

As Kepler continues to monitor transits of hot- Jupiters 
in front of active stars, it will help to better understand 
the stellar cycles. If the Kepler mission is extended, the 
long term photometry will enable it to produce star-spot 
maps and learn more about spot mean lifetimes and pho- 
tospheric differential rotations. In the case of Kepler-17, 
we may be able to measure the complete activity cycle 
for this star and to compare it to another well-know G- 
dwarf: the Sun. 

7.2.3. Impact of stellar variability of Kepler-17 on the 
system parameters 

When the planet transits in front of stellar spots, its 
transit shape deviates from the averaged phase-folded 
lightcurve. The effect of occulted stellar spots on the 
shape of the transit lightcurve is observed in the resid- 
uals from the best-fit transit model of the phase-folded 
light curve (see Figure [2|) . Since the stellar activity in- 
fluences the transit lightcurve profiles, the planetary pa- 
rameters we derive from these profiles are likely to be 
affected. This is a well known proble m for planets tran- 
siting in front of va riable stars (e.g. iCzesla et al.l 120091 : 
iDesert et al.l 1201 lal) . Importantly for the present study, 
the variability affects the stellar density that we assume 
a fixed value for ou r determination of the stellar parame- 
ters (see Section|4]). ICzesla et al.l (2009) propose to fit the 
lower envelope of the transit lightcurve to recover more 
realistic transit parameters. This assumes that dark stel- 
lar structures dominate over bright faculae. In the case 
of Kepler-17b, we cannot exclude the possibility that ev- 
ery transit is affected by dark or bright stellar regions 
so that a priori no individual transit light curve can be 
used as representative of an unaffected profile. Further- 
more, because of the stroboscopic effect described above, 
the phase-folded transit lightcurve possesses combined 
pattern distortions that prevents the use of its lower en- 
velope to derive more accurate parameters. 

We note that the combination of the residuals for the 
eight transit lightcurves presented in Figure [TT] is very 
similar to the total residuals plotted in Figure [5] and ex- 
hibits a symmetrical structure. This symmetrical pattern 
has a shape which is consistent with what would be ex- 
pected from an oblate planet but with an amplitude ten 
to a hundred times larger than expected for a Jupiter 
size planet (Scagcr & Hui 2002; Barnes fc Fortnevll2003l : 
iCarter ^PWinni,2010, ) . This pattern represents a remark- 
able coincidence, in that it seems quite symmetric but it 
is only obtained when all 8 time series are combined, and 
does not appear in the residuals for any single time series. 



Discovery and characterization of Kepler- 17b 



7.3. Atmospheric constraints for Kepler- 17b 

The occuhation depths measured in each bandpass are 
combined and turned into an emergent spectrum for the 
planet. The observed flux of the planet in each band- 
pass corresponds to the sum of the reflected light and 
the thermally emitted light. The combined occultation 
observed in the Kepler bandpass provides a measure of 
the planet's geometric albedo. The occultation depth of 
58 ±8 ppm corresponds to Ag — 0.10 ±0.02. Assuming a 
Lambertian criterion Ab < 1.5 x Ag, and assuming that 
the geometric albedo measured in the Kepler bandpass 
is the same at every wavelength, we infer an upper limit 
to the Bond albedo of Ab < 0.18 at the 1-cr level. This 
result is conservative in the sense that the true albedo 
may be even lower, if some of the occultation signal we 
have measured in the Kepler bandpass is due to thermal 
emission rather than reflected light. This is expected for 
the temperature regime of hot-Jupiters. Therefore, the 
Bond albedo could be well below 0.18. 

We estimate the thermal component of the planet's 
emission in the two Spitzer band-passes from the oc- 
cultation depths measured at 3.6 and 4.5 ^m (see Ta- 
ble [5]). We assume that the planetary emission is well 
reproduced by a black-body spectrum and translate the 
measured depth of the secondary eclipse into brightness 
temperatures. We use t he PHOENIX atmospheric code 
(jHauschildt et all 119990 to produce theoretical stellar 
models for the star Kepler-17. Taking the Spitzer spec- 
tral response function into account, the ratio of areas of 
the star and the planet and the stellar spectra, we derive 
the brightness temperatures that best-fit the observed 
eclipse depths measured in the two IRAC bandpasses. 
The brightness temperature calculated this way result in 
T3.6 ^m=1880 ± 110 K and T4.5 /zm=1770 ± 150 K. 

We compare our data to the h o t- Jup iter atmospheric 
model described in lFortnev et al.l (I2008D which has been 
used for a variety of close-in planets (|Fortnev et al.|[2005l 
1200 61). We aim at broadly distinguishing between differ- 
ent classes of model atmospheres, given that there are no 
prior constraints on the basic composition and structure. 
The atmospheric spectrum calculations are performed for 
ID atmospheric pressure-temperature (P-T) profiles and 
use the equilibrium chemical abunda n ces, at solar metal- 
licity, described in ILoddersI (J2002D : ILodders fc FeglevI 
(poo 6;) . This is a self-consistent treatment of radia- 
tive transfer and chemical equilibrium o f neutral species. 
The o pacity database is described in iFreedman et al.l 
()2008() . The models are calculated for various values of 
the dayside-to-nightside energy redistribution parameter 
(/) and allow for the presence of TiO at high altitude, 
which may play the role of an absorber and which likely 
lead to an inversion of the P-T profile. 

We compare the data to model predictions and select 
models with the best reduced x^ ■ We note that this 
is not a fit involving adjustable parameters. All models 
shown assume redistribution of absorbed flux over the 
dayside only, favoring low values for /. Since our obser- 
vations are only weakly constraining, the 3.6 and 4.5 /ini 
depth ratios are best fitted by a model without a temper- 
ature inversion, although the model with a temperature 
inversion is a reasonable fit as well. In the model with no 
inversion the optical opacity sources are neutral atomic 
Na and K, while in the inverted model the optical opacity 



is dominated by gaseous TiO and VO molecules. The rel- 
atively large occultation depth in the Kepler bandpass is 
not matched by these simple models. This may indicate 
that silicate clouds (as is seen in L-type brown dwarfs, 
but not modeled here) may lead to a higher amount of 
scattered flux from the planet. It may also indicate that 
optical opacity sources are overestimated, such that more 
planetary thermal flux is being seen from deeper layers. 
The growing sample size of Kepler occultation detections 
may show trends with planetary T^s that could help to 
cl arify this i ssue. 

IHartmaiil (|2010D shows that there is a correlation be- 
tween the surface gravity of hot Jupiters and the activity 
levels of the host stars, such that high surface gravity 
planets tend to be found around high-activity stars. The 
position of Kepl er- 17b on such a, diag ram is consistent 
with this trend. 'Knutson e t al.l ()2010D show that there 
could also be a correlation between the host star activity 
level and the thermal inversion of the planetary atmo- 
sphere. In this picture, the strong XEUV irradiation 
from the active stellar host of a hot- Jupiter depletes the 
atmosphere of chemical species responsible for producing 
inversions. The log-R^j^ of Kepler-17 is consistent with 
the case of no thermal inversio n in the framework devel- 
oped by iKnutson et al.l (I2010D. We also evalua ted the 
empirical index defined bv IKnutson et al.l ()2010D . which 
could be correlated with the presence of a thermal inver- 
sion. Using the same definition, we find indices of 0.026% 
above and 0.034% below the predicted black-body fluxes 
in the 3.6 and 4.5 micron bands. These indices suggest 
that the atmosphere of Kepler- 17b is consistent with a 
non-inverted profile, which is in agreement with the re- 
sults of our present study. 
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TABLE 1 

Relative Radial- Velocity Measurements of Kepler- 17 



HJD 


Phase 


RV 


±0"RV 


BS 


io-BS 


(days) 


(cycles) 


(ms-l) 


(ms-i) 


(ms-i) 


(ms-i) 


2455430.839929 


313.013 


-142.61 


67.08 


-77.7 


114.6 


2455431.812259 


313.667 


325.71 


49.17 


-23.8 


72.4 


2455432.819992 


314.345 


-300.32 


90.74 


42.2 


68.9 


2455442.794995 


321.060 


-119.08 


43.22 


45.6 


68.2 


2455443.765027 


325.080 


528.51 


198.75 


483.3 


276.7 


2455448.768015 


327.760 


-158.75 


51.43 


-110.9 


218.2 


2455452.748541 


341.195 


358.36 


61.44 


-50.3 


67.7 


2455472.707059 


346.565 


-413.96 


25.96 


-113.8 


80.9 


2455480.686032 


355.957 


138.23 


56.24 


-116.1 


70.4 


2455494.639064 


357.302 


83.68 


36.83 


-47.9 


75.7 


2455496.636013 


372.077 


-511.48 


60.28 


-36.6 


75.3 


2455518.588037 


372.747 


-205.45 


56.70 


1.2 


67.6 


2455522.579531 


374.764 


417.16 


28.10 


-121.9 


34.2 
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Fig. 1. — Kepler raw long-cadence lightcurve of Kepler-17, from Quarter to Quarter 6. The drop of flux occurring every 1.48 days 
corresponds to the planetary transits. The stellar activity shows 3% flux variation and with a period of 11.89 days. 
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Fig. 2. — Top panel: Kepler transit lightcurve of Kepler- 17b phase-folded, normalized and binned by 30 seconds. The best-fit model of 
this lightcurve is overplotted in red. The residuals from the best-fit trans it ligh tcurve are plotted in the bottom panel. Occulted stellar 
spots produce a symmetrical pattern in phase in the residuals (see Section l7.2.1ll . 
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Fig. 3. — The left panel shows a I-band natural seeing image of Kepler-17 taken with the Lick l-meter telescope (1.2'xl.2'). The right 
panel shows a J-band image of Kepler-17 taken with the wide field camera (WFCAM, I'xl') on the United Kingdom Infrared Telescope 
(UKIRT). These images confirm that there are no companions from 2" to 5" from the star's center, down to a limit of 19th magnitude. 
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Fig. 4. — Orbital solution for Kcpler-17. Top panel: the observed radial velocities obtained with HET HRS are plotted together with the 
velocity curve for a circular orbit with the period and time of transit fixed by the photometric ephemeris. The thick dashed line corresponds 
to the best fit and the thin dashed lines correspond to models with itlo- error on the parameter K. The 7 velocity has been subtracted from 
the relative velocities here and in Table [T] and thus the center-of-mass velocity for the orbital solution is by definition. Bottom panel: 
The velocity residuals from the orbital solution. The rms of the velocity residuals is 52 ms"'^. 
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TABLE 2 
Warm-Spitzer OBSERVATIONS. 



AOR Wavelength (fim) Obs. Date (UT) Select, points Depth (%) Weighted. Avg. depth Bright. T (K) 



1 

3 

2 
4 


40252416 
40251904 
40252160 
40251648 


3.6 
3.6 

4.5 
4.5 


2010-08-29 
2010-09-16 
2010-08-30 
2010-09-18 


2266 
2301 
2315 
2313 


r, oii+U.U4b 

U.Jii_Q QgQ 
U.ZCiO_Q Q3Q 

n 970+0-040 

U.^^U_QQgQ 

•^4R+0060 


0.250 ± 0.030% 
0.310 ± 0.035% 


1880 ±110 
1770 ± 150 
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Fig. 6. — Phase-folded and normalized Kepler occupation lightcurve of Kepler-17b binned by 15 minutes (top panel). The best-fit model 
of the occultation in the Kepler bandpass is overplotted in blue and the residual from this best-fit are shown in the bottom panel. 
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Fig. 7. — Spitzer occultation lightcurves of Kepler-17 observed in the IRAC band-pass at 3.6 (top) and 4.5 fim (bottom). Two visits 
were obtained for each Spitzer bandpass. Top panels : raw and unbinned lightcurves. The red solid lines correspond to the best-fit m odel s 
which include the time and position instrumental decorrelations as well as the model for the planetary occultations (see details in Sect. 15. 2t . 
Bottom panels ; corrected and normalized occultation lightcurves with their best-fit models (in red). The data are binned in 25 minutes 
intervals (50 data points per bin). 
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Fig. 8. — Upper mass limits for a hypothetical perturber as a function of pcrturbcr period normalized to the period of Kepler-17b, Pc/Pt- 
The dark gray zone corresponds to the region of dynamical instability, the white zone corresponds to the region where the presence of a 
companion with a minimum mass would be permitted by the current limits on TTVs. 
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Fig. 9. — Stellar phase folded Kepler normalized long-cadence lightcurves of Kepler- 17 from Quarter 3. The flux is normalized to its 
median value measured over the whole quarter. The stellar rotation period is 8 times the planet's orbital period. This integer period 
ratio allows one to see the transits at the same stellar phase modulo 8. The stellar variability exhibits 3% flux variation with a period of 
approximately 11.9 days. 
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Fig. 10. — Lomb-Scargle periodograms of the Kepler photometry (in black) and of the observed minus computed (O-C) mid-transit times 
(in red) magnified by 30 for comparison purposes. The periodogram from the photometry shows two main peaks at 5.95 and 11.9 days 
while the O-C periodogram exhibits a dominant peak at 6 days. 
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Fig. 11. — Left: Sequence of combined and binned transit lightcurves, with the best-fit model presented in Figure[2]and overplotted in 
red. Each co-added transit corresponds to the combination of 22 individual transits that occurred at epochs modulo eight planetary orbital 
periods. The lightcurves are binned by 100 seconds and they are shifted vertically for display purposes. Each combination of individual 
transits allows us to increase the SNR and to demonstrate that the same spots are occulted during several consecutives transits and epochs. 
The overall combination of these eight transit lightcurves gives the final curves presented in Figure [2] Occulted stellar spots are revealed 
in the combined curves since the stellar rotation period is eight times the planet's orbital period. The same spots are crossed every eight 
transits at a similar orbital phase. Right: Residuals of the best-fit model subtracted from each individual combined lightcurve modulo 
8. The vertical dashed lines correspond to the beginning and to the end of the transits. Five occulted stellar spots are indicated on the 
residuals (A, B, C, D and E) as they appear transits after transits at phase positions expected from the stellar rotation period. This implies 
that the projected spin-orbit angle, A, is very close to for this system. The combination of the residuals of the eight transit lightcurves is 
similar to the total residuals plotted in Figure |2] and exhibits a symmetrical structure. 
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Fig. 12. — Lifetime of stellar spots occulted during transits. The top panel correspond to the complete map constructed from Quarter 
1 to 6 data. The bottom panel is a zoom of the top panel. This time-map image corresponds to the scatter (S/N) measured in a sliding 
box of duration twice the ingress time across the residuals in each transit epoch, relative to mid-transit (Y-axis) as a function of the epoch 
(X-axis). The regions in blue correspond to times when the sections of individual transit lightcurves are deeper than the averaged transit 
lightcurve, whereas the regions in red correspond to less deep sections attributed to cold stellar spots. For display purposes, we used an 
interpolation with an output sampling of 600 samples in the epochal direction and 300 samples in the transit phase direction. 
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Fig. 13. — Day side planet-to-star flux ratios as function of the wavelength for two atmospheric models IFortnev et al.|[2008l ). The black 
filled diamonds and their error bars correspond to the Kepler and Spitzer observations. The dataset is best fitted with models assuming 
full redistribution of the energy to the planetary day-side. The green model represent a non-inverted atmosphere (no TiO). The blue model 
represents an inverted atmosphere (with TiO). The target is slightly better fitted by the non-inverted model. 
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TABLE 3 

System Parameters for Kepler- 17 



Parameter 


Value 


Notes 


Object of interest 






Kcplor Input Catalog number (KIC) 


10619192 




Kepler Object of Interest number (KOI) 


203 




RA (J2000) 


19 : 53 : 34.86 




Dec (J2000) 


-1-47 : 48 : 54.0 




Kepler magnitude 


14.14 




J magnitude 


12.99 




r magnitude 


14.08 




Transit and orbital parameters 






Orbital period P (d) 


1.4857108 ±0.0000002 


A 


Transit Ephemeris (BJDutc) 


2455185.678035inooo26 


A 


Scaled semimajor axis a/R* 


5.48 ±0.02 


A 


Scaled planet radius Rp/R* 


n iQnc(i -1-0.00022 


A 


Impact parameter b = a cos J/R* 


n 9RS+0014 

0.268_(,g^2 


A 


Orbital inclination i (deg) 


87?2±0.15 


A 


Orbital semi-amplitude K (ms^^) 


419.5lli;l 


A,B 


^/ecos(w) 


0.008lO-°ll 


A,B 


v^sm(ui) 


-0-OS^-oml 


A,B 


Orbital eccentricity e 


< 0.011 


A,B 


Transit duration (d) 


n nn/ioc+0. 00007 

0.09485_o 00007 




Observed stellar parameters (fixing the surface 


gravity) 




Effective temperature T^a (K) 


5630 ± 100 


C 


Metallicity [M/H] 


0.3 ±0.1 


C 


Projected rotation -usini (kms^^) 


4.7 ±1.0 


C 


Absolute Systemic Radial Velocity 7 (kms^^) 


-23.82 ±0.10 


C 


Mean radial velocity (kms~l) 


-1-0.40 ±0.10 


B 


Quadratic limb darkening coefficient (ui) 


0.405 ± 0.007 


A 


Quadratic limb darkening coefficient (U2) 


0.262±0:«» 


A 


Derived stellar parameters 






Mass M^,{Mq) 


1.06 ±0.07 


C,D 


Radius R*{Rq) 


1.02 ± 0.03 


C,D 


Surface gravity logg* (cgs) 


4.43 ± 0.02 


A 


Age (Gyr) 


3 ±1.6 


C,D 


Planetary parameters 






Mass Mp (Mj) 


2.45 ±0.11 


A,B,C,D 


Radius Rp (Rj, equatorial) 


1.31 ±0.02 


A,B,C,D 


Density pp (gcm~3) 


1.35 ±0.08 


A,B,C,D 


Surface gravity loggp (cgs) 


3.55 ± 0.02 


A,B 


Orbital semimajor axis a (AU) 


09591+000037 
U.U^Oyi^o. 00036 


E 


Equilibrium temperature Toq (K) 


1570 ± 200 


F 


Geometric albedo Ag 


0.10 ±0.02 


A 



Note. — 
A: Based on the photometry. 
B: Based on the radial velocities. 
C: Based on an analysis of the TRES spectrum. 
D: Based on the Yale-Yonsei stellar evolution tracks. 
E: Based on Newton's version of Kepler's Third Law and total mass. 
F: Assumes Bond albedo — 0. and complete redistribution. 
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